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Abstract 

The BaBar experiment and the PEP-II accelerator at SLAC started to take data on May 26, 
1999. By the time of this conference, the recorded integrated luminosity was 20 fb _1 , of which 8 
fb _1 were analyzed to provide a first set of physics results. This talk reviews the first measurement 
of sin 2/3 and the study of B meson decays to charmonium modes and 2-body charmless decays. 
Complementary results were presented by the BaBar Collaboration at this conference and are 
reviewed in |fll|. 



Contribued to the Proceedings of the International Conference On CP Violation Physics, 

18-22 Sep 2000, Ferrara, Italy 



Stanford Linear Accelerator Center, Stanford University, Stanford, CA 94309 



Work supported in part by Department of Energy contract DE-AC03-76SF00515. 



1 Physics at BABAR 



The BaBar experiment |Q] is designed in order to perform a comprehensive study of CP violation 
in the B system. The main goal is to investigate if CP violation can be fully explained within the 
Standard Model by the imaginary phase of the Cabibbo-Kobayashi-Maskawa matrix, or if other 
sources of CP violation should be looked for elsewhere Q. 

The strategy is to over-constrain the Unitarity Triangle by measuring at the same time its sides 
and angles. 

The sides can be extracted from measurements of |V^V^| in the study of B°B° mixing, V u b 
from the study of decays of B mesons in charmless final states, and V c f, from studies of the decay 
B -> D*lv. 

The determination of the angles requires to measure time dependent asymmetries for the decay 
of a neutral B meson into a CP eigenstate, fcp, of the form 

A r(B (t)^fcp)-T(W(t)^f CP ) 
CP{ r(B°(t) -> f CP ) + T(B°(t) - fcp) ' 

For B° decays with only one diagram contributing to the final state 

A C p{t) = -ImAsin(AmAt), (2) 

where A = p/qA(f)/A(f), A(f) and A(f) are the decay amplitudes and p and q are the complex 
coefficients that relate mass and flavour eigenstates of the B mesons ||. For some decay modes, 
ImA is directly and simply related to the angles of the Unitarity Triangle. Of particular interest 
are the decays of the neutral B mesons in final states containing a charmonium state and a neutral 
kaon. For these decays 

Acp(t) = ± sin(2/3) sin(AmAt). (3) 
Experimentally, the measurement of the CP asymmetry Acp requires three basic steps: 

• the exclusive reconstruction of one B meson (Bcp) in a CP eigenstate (e.g. B° — > J/ipK s ); 

• the measurement of the time At elapsed between the decay of the two neutral B mesons 
produced in the T(45) decay; 

• the tagging of the CP eigenstate at the production time, inferred from the tagging of the 
other B meson (B tag ) produced coherently in the T(45) decay. 



2 The experimental apparatus 

The PEP-II storage ring and the BaBar detector have been designed specifically to measure time- 
dependent CP asymmetries in the neutral B meson decays in the most effective way. 



2.1 The PEP-II accelerator 

PEP-II is an asymmetric B factory consisting of two concentric rings: a high energy ring (HER) 
with a 9.0 GeV electron beam and a low energy ring (LER) with a 3.1 GeV positron beam. The 
energy of the two beams has been chosen in order to achieve a center of mass energy of 10.58 GeV, 
the mass of the T(4S) resonance, and a boost of 0.56 of the T(4S) along the z axis. The boost 



Table 1: Typical reconstruction efficiencies and pion misidentification probabilities for electrons, 
muons and kaons in BaBar. 



Particle 


Efficiency 


7T contamination 


electrons 


92% 


0.3% 


muons 


75% 


2.5% 


kaons 


85% 


5% 



corresponds to an average separation of 260 fim between the decay vertices of the two B mesons 
produced in the T(4S) decay. The design luminosity of the machine is 3x 10 33 s _1 cm -2 , which 
corresponds to about 30 million BB pairs produced per year. 

PEP-II started its activity in factory mode on May 26, 1999 and reached outstanding per- 
formance very rapidly. By the time of this conference the achieved peak luminosity was 2.56 
xl0 33 s _1 cm -2 , which corresponds to 85% of the design value. The design daily recorded luminos- 
ity of 135 pb _1 /day was routinely achieved and exceeded starting in June 2000, with a record of 
151 pb _1 accumulated in one day. By the time of this conference the total integrated luminosity 
delivered by PEP-II was 19.5 fb _1 , of which 17.5 recorded at the T(AS) resonance and 2 fb _1 at 
a center of mass about 40 MeV below resonance. Thanks to a data taking efficiency higher than 
95%, BaBar accumulated more than 18 fb" 1 on tape. 

2.2 The BABAR detector 

BaBar is a cylindrical detector built around the interaction point of PEP-II. It consists of several 
concentric subdetectors. From the inside out there are: 

• a silicon vertex tracker, consisting of 5 layers of double-sided silicon strips; 

• a cylindrical drift chamber, filled with a low density gas mixture (Helium-Isobutane: 80%- 
20%) to minimize multiple scattering; 

• the DIRC detector, the main particle identification device of the experiment, which consists 
of 144 quartz bars; 

• an electromagnetic calorimeter, consisting of 6580 Csl crystals; 

• a 1.5 Tesla superconducting magnet; 

• an instrumented flux return (IFR), consisting of 19 layers of RPCs, for muon identification 
and Kl reconstruction. 

The combined tracking resolution of the BaBar experiment can be parameterized as 

( = (0.0015p T ) 2 + 0.005 2 . (4) 

V PT J 

The performances of electron, muon and kaon identification are summarized in Table |l]. Of 
particular interest is the good kaon-pion separation capability of the detector obtained from the 
measurement of the Cherenkov angle by the DIRC detector. The separation is better than 3a up 
to momenta of 3.5 GeV/c as illustrated in Figure |l|. 

A more detailed description of detector design and performance can be found in Q. 
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Figure 1: (a) The Cherenkov angle and (b) K-ir separation as functions of momentum for single 
tracks in a D° control sample. The separation is an average over all polar angles. 



2.3 The data set 

Unless otherwise specified, all the results presented in this talk are based on a data set of 7.7 fb _1 
recorded at the T(AS) resonance and 1.2 fb _1 recorded off peak. These data were recorded by the 
BaBar detector between January and June 2000. 

Since the analyses are still in progress, all results are preliminary. 



3 Study of B — > charmonium decays 

A good understanding of B meson decays into final states containing a charmonium resonance is 
a prerequisite to an analysis of CP violation in the B system. In this section, an extensive set of 
measurements of inclusive and exclusive branching fractions of B mesons into charmonium decays 
is briefly summarized. More details can be found in || and ||. 

3.1 Inclusive analysis 

J/ip candidates are formed from pairs of oppositely charged tracks both identified as leptons. 
Figure H shows the invariant mass distributions for J/ip — * e + e~ and J/ifi — * candidates after 

continuum subtraction. The number of J/ip mesons is extracted from the mass distribution using 
a probability density function (p.d.f.) derived from a simulation that includes final state radiation 
and bremsstrahlung. A total of 4920 ± 100 ± 180 J/tp -> e+e" and 5490 ± 90 ± 90 J/ip /z+ \T 
events were reconstructed. 

The 4>(2S) was reconstructed in its decays to e + e~, [J, + fi~ and J/V'vr + 7r~. For the modes fi~ 
and e + e~, the numbers of signal events are extracted from the mass distribution using a p.d.f. 
from simulation. A total of 131 ± 29 ip(2S) — > e + e~ and 125 ± 19 ip(2S) — > /x + /U~ events were 
found. For the decays J/^7r + 7r _ , the yields are extracted from a fit to the distribution of the mass 
difference between the ip{2S) and the J/ip in order to reduce the impact of the radiative tail and 




Figure 2: Mass distribution of J/tp candidates to (a) electron pairs and (b) muon pairs after the 
continuum subtraction. 

the mass resolution of the J/tp candidate. A total of 126 ±44 events were found in candidates with 
J/tp — > e + e~ , and 162 ± 23 events for J/tp —* 

Xc mesons are reconstructed through the radiative decay Xc — ► jJ/ip ■ The selected photon 
candidates must have an energy between 0.20 and 0.65 GeV in the CM system and be isolated from 
nearby hadronic showers. The substantial combinatorial background from ir° decays is reduced 
by rejecting any photon that, when combined with any other photon, gives an invariant mass 
compatible with a ir°. 

As for the ip(2S), the numbers of events are extracted by fitting the distribution of the mass 
difference between the Xc and J/tp candidates (Figure ||). We simultaneously fit for a Xci and a 
possible Xc2 component, assuming the resolution of the two resonances to be the same and fixing 
the difference between the Xc2 and the Xci masses to the Particle Data Group value. We found 
129 ± 26 ± 13 Xci and 3 ± 21 Xc2 events in which J/tp decayed into e + e~ , and 204 ± 47 ± 12 Xci 
and 47 ± 21 Xc2 events for the fi + pT mode. 

Due to the preliminary status of the analysis, we chose not to quote a branching ratio for B — > 
J/tp X at this time. The yields obtained for the decay B — > J/tp X are used, together with the 
already precise value of the branching fraction for this decay reported in Q, as a normalization in 
the calculation of the branching ratios of B — > tp(2s)X and B — > XcX. The results, summarized in 
Table g, are very competitive compared with the world averages. 

3.2 Exclusive analysis 

We considered the exclusive B decay channels listed in Table ||. 

In these analyses, the reconstruction of the charmonium decays is very similar to the one used 
for inclusive analysis, but with looser requirements for the lepton identification. In the case of the 
decay J/tp — ► e + e~ , we apply a procedure to add photons which are close to the electron tracks in 
order to reduce the impact of bremsstrahlung on the reconstruction efficiency. 

K® — > 7r + 7r~ candidates are formed from pairs of oppositely charged tracks which have an 
invariant mass between 0.489 and 0.507 GeV/ c 2 . In addition we require that the K® be consistent 




Figure 3: Mass difference between the \ c and J/ip candidates, where \c — ► jJ/ip an d (a) J/ip — > 
e + e~ and (b) J/ip — ► /J- + fi~. 



Table 2: Summary of inclusive branching ratios measured with respect to B — ► J/ipX. The 
results are combined and multiplied by the PDG value for ,6(-B — > J/ipX) to obtain inclusive 
branching fractions (%). 





Branching ratio relative to J/ipX 


Branching fraction (%) 




e + e 

Stat Sys 


Stat Sys 


Common 

Sys (%) 


Combined 

Stat Sys 


Xci 7«W 


0.28 ± 0.06 ±0.03 


0.38 ±0.05 ±0.02 


±7.3 


0.39 ±0.04 ±0.04 


Xc2 90% CL 


< 0.16 


< 0.28 


±8.1 


< 0.24 


^(25) -»• e+e- 

ip(2S) -► ^7T J/0 
Combined Y>(2S) 


0.22 ±0.05 ±0.01 
0.18 ±0.06 ±0.01 


0.23 ±0.04 ±0.01 
0.22 ±0.03 ±0.00 


±13.6 
±9.9 


0.26 ±0.03 ±0.04 
0.24 ±0.03 ±0.03 
0.25 ±0.02 ±0.02 



with having originated from the J/ip vertex. 

K® — ► 7r°7r° candidates are required to have a mass between 0.470 and 0.525 GeV/c 2 and an 
energy greater than 0.8 GeV. A ir° decay to two photons is observed in the EMC either as a single 
neutral cluster with substructure or as two distinct 7 clusters. The most probable decay point of 
the Kg is determined after refitting the two ir° mesons at several points along the path defined by 
their summed momentum vector and the J/tp vertex. 

We reconstruct K*° decays to K + n~ and K^n , and K* + decays to Kg-K + and K + tt°. In all 
cases the candidate K* is required to have an invariant mass within 0.075 GeV/ c 2 of the nominal 
value. 

Charged kaons are identified by the Cherenkov angle measured by the DIRC detector and/or 
by the specific ionization measurement (dE/dx) in the DCH, depending on particle momentum. 

To isolate the signal for each mode we use the variables AE, the difference between the re- 
constructed and expected B meson energy measured in the center-of-mass frame, and toes> the 
beam-energy substituted mass. These variables are defined as: 

™es = ^jEf-pf, (5) 
AE = E* B -E* b , (6) 

where E£ is the beam energy in the center-of-mass, and E* B and p* B are the energy and momentum 
of the reconstructed B meson in the center-of-mass. 

When deriving branching fractions we have used the secondary branching fractions and their 
associated errors published by the Particle Data Group M . 

We determine the number of BB events from the difference in the multi-hadron rates on and off 
the T(4S) resonance, normalized to the respective luminosities. The efficiencies for each mode have 
been obtained from Monte Carlo simulations complemented with measurements of tracking and 
particle identification efficiencies extracted from data. The shape of the beam-energy substituted 
mass distribution is parameterized for each mode with the sum of a Gaussian and the ARGUS 
function M. 

Figure^ shows the mES an d AE distributions of the candidates. In Table || we present the 
yields and measured branching fractions for the individual exclusive modes. Figure [| shows the 
measured branching fractions compared to the values compiled by the Particle Data Group M . 



Table 3: Exclusive B meson decay modes. 



Channel 


Secondary decay mode(s) 
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Figure 4: Distributions of candidate events in mES an d AE. The upper plot shows the B modes 
and the lower plot the B + modes. 



Table 4: The yields and measured branching fractions for exclusive decays of B mesons involving 
charmonium. The errors on the yields are only statistical. For the branching fractions, the first 
error is statistical and the second systematic. 



Channel Yield Br. Frac/10 



IT 



K° 7r + 7r~ 93 ± 10 10.2 ±1.1 ±1.3 

K° -> 7r°vr 14 ± 4 7.5 ± 2.0 ± 1.2 

B + ^J/i)K+ 445 ± 21 11.2 ±0.5 ±1.1 

B°^J/^K*° 188 ± 14 13.8 ± 1.1 ± 1.8 

B + ^J/^K* + 126 ± 12 13.2 ±1.4 ±2.1 

B°^ifj(2S)K° 23 ± 5 8.8 ±1.9 ±1.8 

B + ^(2S)K + 73 ± 8 6.3 ±0.7 ±1.2 

B+^XdK+ 44 ± 9 7.7 ±1.6 ±0.9 
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Figure 5: Summary of branching fraction measurements for charmonium + K channels and 
comparisons with the PDG 2000 values. 



4 Measurement of sin 2/3 

The measurement of sin2/3 requires five basic steps: exclusive reconstruction of the CP mode; 
measurement of time elapsed between the decays of the two B mesons; tagging of the flavour of the 
B meson at production time; measurement of the dilution factors from the data; extraction of the 
value of sin2/3 from a fit to the amplitude of the CP asymmetry. The following subsections briefly 
describe the analysis strategy and the result. For a more detailed description see || and references 
therein. 

4.1 Reconstruciton of the CP mode 

For a first measurement, we targeted only decays characterized by relatively high accessible branch- 
ing fractions (~10 _4 -10 -5 ), low backgrounds (<10%), and reliable theoretical prediction (i.e., no 
penguin pollution). Only the decays B — > J/ipK s and B — * ip(2S)K s have been used so far. The 
procedure used for the exclusive reconstruction of the B meson is described in [||. 

The yields and purities obtained in each channel analysing 9 fb _1 are reported in Table [|. The 
modes B — ► XcK s , B — ► J/tpK^ and B — > J/ipK*° are still under investigation and will be included 
in the analysis in the future. 

4.2 Time resolution function 

The resolution of the At measurement is dominated by the z resolution of the tagging vertex. The 
three-momentum of the tagging B and its associated error matrix are derived from the fully recon- 
structed Bqp candidate three momentum, decay vertex and error matrix, and from the knowledge 
of the average position of the interaction point and the T(4S) four-momentum. This derived Bt ag 
three-momentum is fit to a common vertex with the remaining tracks in the event (excluding those 
from Bqp). In order to reduce the bias due to long-lived particles, all reconstructed V° candidates 
are used as input to the fit in place of their daughters. Any track whose contribution to the x 2 
is greater than 6 is removed from the fit. This procedure is iterated until there are no tracks 
contributing more than 6 to the x 2 or until all tracks are removed. 

The time resolution function is described by the sum of two Gaussian distributions: 



Monte Carlo simulation indicates that 70% of the events are in the core Gaussian, which has a 
width G\ ~ 0.6 ps. The wide Gaussian has a width oi ~ 1.8 ps. Tracks from forward-going charm 
decays included in the reconstruction of the Btag vertex introduce a small bias, b\ ~ — 0.2ps, for 
the core Gaussian. 

Table 5: The CP sample used for the sin2/3 measurement: yields and purities obtained using 9 




(7) 



fb- 1 . 



Decay mode 



Yield Purity 



J/^K S {K S ■ 
J/^K S (K S ■ 
^(2s)K s (K t 



TT+TT") 124 ± 12 96% 
A ) 18 ± 4 91% 
-> vr+vr-) 27 ± 6 93% 



Table 6: Parameters of the resolution function determined from the sample of events with fully- 
reconstructed hadronic B° candidates. 



Parameter 


Value 


Si (ps) 


-0.20 ±0.06 


from fit 


Si 


1.33 ±0.14 


from fit 


U (%) 


1.6 ±0.6 


from fit 


fl (%) 


75 


fixed 


S 2 (ps) 





fixed 


s 2 


2.1 


fixed 



A small fraction of events (~1%) have very large values of Az, mostly due to vertex reconstruc- 
tion problems. This is accounted for in the parameterization of the time resolution function with 
a very wide unbiased Gaussian with fixed width of 8ps. 

In likelihood fits, we use a parameterization of the error a At calculated from the fits to the two 
B vertices for each individual event. We introduce two scale factors S± and S 2 for the width of the 
narrow and the wide Gaussian distributions (a± = Si X a At and a 2 = S 2 x a At) to account for the 
fact that the uncertainty on At is underestimated due to effects such as the inclusion of particles 
from D decays and possible underestimation of the amount of material traversed by the particles. 
The scale factor Si and the bias <5i of the narrow Gaussian are free parameters in the fit. The scale 
factor S 2 and the fraction of events in the wide Gaussian, f 2 , are fixed to the values estimated from 
Monte Carlo simulation by a fit to the pull distribution (52 = 2.1 and f 2 = 0.25). The bias of the 
wide Gaussian, 5 2 , is fixed at Ops. The last free parameter in the fit is f w , the fraction of the very 
wide Gaussian that describes the tail of the resolution. 

Because the time resolution is dominated by the precision of the Bt ag vertex position, we find 
no significant differences in the Monte Carlo simulation of the resolution function parameters for 
the various fully reconstructed decay modes. This allows us to determine the resolution function 
parameters with the relatively high-statistics fully-reconstructed B° data samples. The differences 
in the resolution function parameters in the different tagging categories (see below) are also small. 

Table || presents the values of the At resolution parameters obtained from a maximum likelihood 
fit to a sample of fully reconstructed B° decays. Further details on the procedure and the results 
can be found in |l(| . The vertex parameters are fixed to these values in the final unbinned maximum 
likelihood fit for sin2/? in the low-statistics CP event sample. 

4.3 B flavor tagging 

Each event with a CP candidate is assigned a B° or B° tag if the rest of the event (i.e., with the 
daughter tracks of the Bqp removed) satisfies the criteria for one of several tagging categories. The 
figure of merit for each tagging category is the effective tagging efficiency Qi = e, (1 — 2wi) 2 , where 
£j is the fraction of events assigned to category i and Wi is the probability of misclassifying the tag 
as a B° or B° for that category. Wi are called mistag fractions. The statistical error on sin2/3 is 
proportional to where Q = J2iQi- 

The algorithm used in the analysis categorizes the events in four different classes: 

Lepton when the flavor of the B meson is identified by the charge of the high momentum lepton 
coming from a semileptonic decay; 



Table 7: Categories of tagged events in the CP sample. The Lepton category is split into Electron 
and Muon categories. 



Tagging Category 


J/^K° S 


^(25)^ 


CP sample 






-TP") 




— ► 7T 


°7T ) 






-7T-) 


(tagged) 


B° 


B° 


all 


B° 


B° 


all 


B° 


B° 


all 


B° 


B° all 


Electron 


1 


3 


4 


1 





1 


1 


2 


3 


3 


5 8 


Muon 


1 


3 


4 











2 





2 


3 


3 6 


Kaon 


29 


18 


47 


2 


2 


4 


5 


7 


12 


36 


27 63 


NT1 


9 


2 


11 


1 





1 


2 





2 


12 


2 14 


NT2 


10 


9 


19 


3 


3 


6 


3 


1 


4 


16 


13 29 


Total 


50 


35 


85 


7 


5 


12 


13 


10 


23 


70 


50 120 



Kaon when the flavor of the B meson is identified by the charge of the kaon coming from the 
hadronization of the s quark coming from the decay b — ► c — ► s; 

NT1 the tag is chosen according to the output of a neural network that exploits other (correlated) 
information, such as secondary lepton charge, slow pions from D* decays, jet charge; 

NT2 similar to NT1 but characterized by a smaller tagging power. 

When an event satisfied the selection criteria for more than one tagging category, it was assigned 
to the tag with the highest discrimination power. 



4.3.1 Measurement of mistag fractions 

The mistag fractions are measured directly in events in which one B° candidate, called the B rec , 
is fully reconstructed in a flavor eigenstate mode. The flavor-tagging algorithms described in the 
previous section are applied to the rest of the event, which constitutes the potential Bf ag . 

Considering the B°B° system as a whole, one can classify the tagged events as mixed or unmixed 
depending on whether the Bt ag is tagged with the same flavor as the B rec or with the opposite 
flavor. The observed fraction of mixed events at the time t is expressed as 

X(At) =w + D/2(l -cos(AmAi)), (8) 

were D = 1 — 2w is the so called dilution f actor . The mixing probability is smallest for small 
values of At = t rec — ttag so that the apparent rate of mixed events near At = is governed 
by the mistag probability. In addition to improving sensitivity to the mistag fraction, this time- 
dependent measurement technique also helps discriminate against backgrounds with different time- 
dependence. 

The extraction of the mistag probabilities is complicated in reality by the possible presence 
of mode-dependent backgrounds. We deal with these by adding specific terms in the likelihood 
functions describing the different types of backgrounds (zero lifetime, non-zero lifetime without 



mixing, non-zero lifetime with mixing). Details are described in ||lj] and |1C]. 

The mistag fractions and the tagging efficiencies obtained by combining the results from max- 
imum likelihood fits to the time distributions in the B° hadronic and semileptonic samples are 
summarized in Table ||. We find a tagging efficiency of (76.7 ± 0.5)% (statistical error only). The 



Table 8: Mistag fractions measured from a maximum-likelihood fit to the time distribution for 
the fully-reconstructed B° sample. The uncertainties on e and Q are statistical only. 



Category 


e{%) 


w (%) 


Q{%) 


Lepton 


11.2 ±0.5 


9.6 ± 1.7 ± 1.3 


7.3 ±0.7 


Kaon 


36.7 ±0.9 


19.7 ± 1.3 ± 1.1 


13.5 ± 1.2 


NT1 


11.7 ±0.5 


16.7 ±2.2 ±2.0 


5.2 ±0.7 


NT2 


16.6 ±0.6 


33.1 ±2.1 ±2.1 


1.9 ±0.5 


all 


76.7 ±0.5 




27.9 ±1.6 
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Figure 6: Variation of the log likelihood as a function of sin2/3. The two horizontal dashed lines 
indicate changes in the log likelihood corresponding to one and two statistical standard deviations. 



lepton categories have the lowest mistag fractions, but also have a low efficiency. The Kaon cat- 
egory, despite having a larger mistag fraction (19.7%), has a higher effective tagging efficiency; 
one-third of events are assigned to this category. Altogether, lepton and kaon categories have an 
effective tagging efficiency Q ~ 20.8%. The neural network categories increase the effective tagging 
efficiency by ~ 7% to an overall Q = (27.9 ± 0.5)% (statistical error only). 

Of the 168 CP candidates, 120 are tagged: 70 as B° and 50 as B°. The number of tagged 
events per category is given in Table [7]. 

4.4 Results 

The maximum-likelihood fit for sin2/3, using the full tagged sample of B°/B° — ► J/ipK® and 
B°/B° -> ip(2S)K% events, gives: 

sin2/3 = 0.12 ± 0.37(stat) ± 0.09(syst) . (9) 

For this result, the B° lifetime and Am^ are fixed to the current best values and At resolution 
parameters and the mistag rates are fixed to the values obtained from data as summarized in 
Tables [6| and ||. The log likelihood is shown as a function of sin2/3 in Figure ^. 

The results of the fit for each type of CP sample and for each tagging category are given in 
Table ||| The contributions to the systematic uncertainty are summarized in Table 10, 



Table 9: Result of fitting for CP asymmetries in the entire CP sample and in various subsamples. 



sample 


sin2/3 


CP sample 


0.12±0.37 


Jli)K% (K" -► TT+vr") events 


-0.10 ±0.42 


other CP events 


0.87 ±0.81 


Lepton 


1.6 ± 1.0 


Kaon 


0.14 ±0.47 


NT1 


-0.59 ±0.87 


NT2 


-0.96 ±1.30 



Table 10: Summary of systematic uncertainties. The different contributions to the systematic 
error are added in quadrature. 



Source of uncertainty 


Error on sin2/3 


r U 
T B 


0.002 


Am d 


0.015 


Az resolution 


0.019 


time-resolution bias 


0.047 


measurement of mistag fractions 


0.053 


different mistag fractions 




for CP and non-CP samples 


0.050 


for B° and B° 


0.005 


background 


0.015 


total systematic error 


0.091 



Table 11: Results of fitting for apparent CP asymmetries in various charged or neutral flavor- 
eigenstate B samples. 



Sample 


Apparent CP-asymmetry 


Hadronic charged B decays 


0.03 ±0.07 


Hadronic neutral B decays 


-0.01 ± 0.08 


J/ipK+ 


0.13 ±0.14 


J/4>K*° (K*° -► K+tt-) 


0.49 ± 0.26 



We estimate the probability of obtaining the observed value of the statistical uncertainty, 0.37, 
on our measurement of sin2/3 by generating a large number of toy Monte Carlo experiments with 
the same number of tagged CP events, and distributed in the same tagging categories, as in the CP 
sample in the data. We find that the errors are distributed around 0.32 with a standard deviation 
of 0.03, and that the probability of obtaining a value of the statistical error larger than the one we 
observe is 5%. Based on a large number of full Monte Carlo simulated experiments with the same 
number of events as our data sample, we estimate that the probability of finding a lower value of 
the likelihood than our observed value is 20%. 

To validate the analysis we use the charmonium control sample, composed of B + — ► J/ipK + 
events and events with self-tagged J/ipK*° (K*° — > K + ir~) neutral B's. We also use the event 
samples with fully-reconstructed candidates in charged or neutral hadronic modes. These samples 
should exhibit no time-dependent asymmetry. In order to investigate this experimentally, we define 
an "apparent CP asymmetry" , analogous to sva2(5 in Eq. |l|, which we extract from the data using 
an identical maximum-likelihood procedure. 

The events in the control samples are flavor eigenstates and not CP eigenstates. They are used 
for testing the fitting procedure with the same tagging algorithm as for the CP sample and, in the 
case of the B + modes, with self-tagging based on their charge. We also perform the fits for B° and 
B° (or B + and B~~) events separately to study possible flavor-dependent systematic effects. For 
the charged B modes, we use mistag fractions measured from the sample of hadronic charged B 
decays. 

In all fits, including the fits to charged samples, we fix the lifetime t b o and the oscillation 
frequency Am^ to the PDG values [jjj . The results of a series of validation checks on the control 



samples are summarized in Table 11. The two high-statistics samples and the J/ipK + sample give 
an apparent CP asymmetry consistent with zero. The 1.9a asymmetry in the J/tpK*° is interpreted 
as a statistical fluctuation. 

Other BABAR time-dependent analyses presented at this Conference demonstrate the validity of 
the novel technique developed for use at an asymmetric B Factory. In particular, the measurement 



of the B°-B° oscillation frequency described in [10] uses the same time resolution function and 



tagging algorithm as the CP analysis, and the B lifetime measurement described in [11 1 uses the 
same inclusive vertex reconstruction technique as the CP analysis. Both result are consistent with 
the world average 0. 

4.5 Current knowledge of the angle (3 

Figure [?] shows the (p, ff) plane, with BABAR's measured central value of sin2/3 shown as two 
straight lines. There is a two- fold ambiguity in deriving a value of (5 from a measurement of 
sin2/3. Both choices are shown with cross-hatched regions corresponding to la and 2a experimental 
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Figure 7: Present constraints on the position of the apex of the Unitarity Triangle in the (p, fj) 
plane. Our result sin2/3 = 0.12 ± 0.37(stat) is represented by cross-hatched regions corresponding 
to one and two statistical standard deviations. The fitting procedure is described in Ref. (T^j. 



uncertainty. The ellipses correspond to the regions allowed by all other measurements that constrain 
the Unitarity Triangle for a variety of representative theoretical parameters. This procedure is 
discussed in detail in ifLq ]. The following set of measurements and model-dependent parameters 
are used: [V^l = 0.0402 ± 0.017; \V ub /V cb \ = (\V ub /V cb \) ± 0.0079; Am d = 0.472 ± 0.017 Tips" 1 ; 
| ex | = (2.271 ± 0.017) x 10 -3 ; the set of amplitudes corresponding to a 95%CL limit of 14.6 Tips -1 
for Ara s ; (\V ub /V cb \) in [0.070, 0.100]; B K in [0.720, 0.980]; f Bd y/B^ in [185, 255] MeV; and £ s 
in [1.07, 1.21]. 

In Figure ||, a summary of all existing direct measurements of sin 2/3 is reported. The PDG 2000 
world average, obtained from the OPAL jl3| and CDF |14| results, is compared with a new average 



that includes the recent ALEPH [fLq |, Belle [16| and BaBar j|] results. The new measurements 
from the two B factories improved the knowledge of sin 2(3 by a factor 2. 



5 Study of B — > charmless 2-body decays 
5.1 Interest of the measurement 

The branching fractions for the rare charmless decays B° — » h + h~ (h = 7T,K) provide important 
information in the study of CP violation. If the tree level diagram dominates the decay, the 7r + 7r~ 
decay mode can be used to extract the angle a of the Unitarity Triangle. If the penguin contribution 



is not negligible, large theoretical uncertainties [17] complicate the analysis and the extraction of a 



will require a full isospin analysis [18]. On the other hand, large penguin amplitude improves the 
prospects for observing direct CP violation as an asymmetry in the decay rates for B° — ► K + ir~ 
and B° —> K~ir + . 

The decay B° — ► K + tt~ is dominated by the b — ► sg penguin amplitude, and provides an 
estimate of the scale of penguin pollution in the 7r + 7r~ decay. Precise measurement of the decay 
rates for tttt and Kir decays is therefore of central importance. 



Direct measurements of sin2(3 



OPAL 



CDF 99 



ALEPH 



Belle 



BaBar 
PDG 2000 
All measurements 



sin2|3 



0.79 ±0.39 ±0.16 



0.84 +0 , 8 „ 4 ±0.16 



+0.43 +0.07 
-0.44 -0.09 

0.1 2 ±0.37 ±0.09 

0.9 + 0.4 
0.49 + 0.23 

f*=Preliminarv) 



Figure 8: Existing direct measurements of sin 2(3. The input measurements are in Refs. |U|], [14|, 

PL m, a. 



5.2 Analysis strategy 

The first step of this analysis consists in selecting decays B — >h + h~ , where /i ± is either a pion 
or a kaon. B decays to charm mesons are not a significant background to charmless two-body 
decays because of the relatively small CM momenta of decay products produced from the b — > c 
transition. The background is dominated by the continuum production of light quarks, e + e~ — > qq 
{q = u,d,s,c), which typically exhibits a two-jet structure in the CM frame. The topology is 
markedly different from that of real B — > h + h~ events, which is more spherically symmetric. This 
difference is exploited using the angle 9s between the sphericity axes, evaluated in the CM frame, 
of the B candidate and the remaining charged and neutral particles in the event. 

Further separation between signal and continuum background is provided by a Fisher discrim- 
inant technique |p!9| . The Fisher discriminant T is calculated from a linear combination of 9 
discriminating variables constructed from the scalar sum of the momenta of all charged and neu- 
tral particles (excluding the candidate daughter tracks) flowing into 9 concentric cones centered 
on the -B-candidate thrust axis in the CM frame. More energy will be found in the cones nearer 
the candidate thrust axis in jet-like continuum background events than in the more isotropic BB 
events. 

Once the decay B — > h + h~ is identified, the separation between different channels is obtained 
using two different approaches: 

• likelihood analysis, in which the K-tt separation is determined on a statistical basis; 

• cut based analysis, in which the K-tt separation is performed track-by-track. 

A good understanding of the particle identification performances is essential in both analyses. 



Table 12: Summary of B — * two-body branching fraction results. Shown are the central fit values 
Ns, the statistical significance, and the measured branching fractions B. For the KK mode, the 90% 
confidence level upper limits are given. The first errors are statistical and the second systematic. 



Mode N s Stat. Sig. ff(lO^) 

~^+^ 29ff^ hTo 9.3^;^i : 

K+7T- 38±1±1 6.7a 12.5±|:g±i 

K+K- 7+1 (< 15) 2.1a < 6.6 




Figure 9: The central value (filled circle) for B(B° — > it + it~) and B(B° — ► K + 7r~) along with the 
no" statistical contour curves for the global likelihood fit, where racr corresponds to a change of n 2 
in — 21og£. 



5.3 Results 

We determine branching fractions for 7r + 7r~ and K + tt~ decays and an upper limit for the K + K~ 



decay using the results of the global likelihood fit. The results are summarized in Table 12. For the 
KK mode we calculate the 90% confidence level upper limit yield and decrease the efficiency by the 
total systematic error (24%) before calculating the upper limit branching fraction. The statistical 
significance of a given signal yield is determined by setting the yield to zero and maximizing the 
likelihood with respect to the remaining variables. 

Figure |9| shows the no likelihood contour curves, where a represents the statistical uncertainty 
only. The curves are computed by maximizing the likelihood with respect to the remaining variables 
in the fit. 



6 Summary 



PEP-II and BaBar had a beautiful start. In their first year of data taking, they have reached 
about 85% of the design peak luminosity, exceeded the design daily luminosity and accumulated 
almost 20 fb^ 1 on tape. 

The first physics results, based on about 8 fb _1 , have been produced. This talk reviews the first 
results obtained by BaBar in the measurement of sin 2/3 and in the studies of B —* charmonium 
and B — ► two body charmless decays. 

The results for the data collected in the entire first run of data taking (about 24 fb _1 ) is expected 
to be available in Spring 2001. 
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